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ABSTRACT 


An attempt heis been made to predict the breakthrough curves for sulphur dioxide in 
a packed bed of solid reactants. An integrated software package has been developed 
to solve the involved differential equation and obtain the breakthrough curves using 
starbase graphics. Ba.sic parameters required for predicting the breakthrough curves 
have been taken from the literature. Breakthrough curves for sulphur dioxide have 
been predicted using two solid reactants, namely, active sodium carbonate and calcined 
magnesium carbonate. For active sodium carbonate, breakthrough curves have been 
predicted for temperatures 150,200,250,300, and 350 "C with an inlet sulphur dioxide 
concentration of 3500 ppm and for different inlet sulphur dioxide concentrations of 
1500, 3500, 6000, and 18000 ppm keeping temperature constant at 200 °C . For calcined 
magnesium carbonate these have been predicted for temperatures 650, 750, and 850 
°C with an inlet sulphur dioxide concentration of 9000 ppm and for concentrations 
2000, 9000, and 22000 ppm keeping temperature constant at 750 °C. Model predicted 
breakthrough curves have been found to match fairly well with the experimental ones, 
except, in the Ccise of those at higher temperatures (350 °C for active sodium carbonate 
and 850 °C for calcined magnesium carbonate). Departure may be attributed to the 
phenomenon of sintering at these temperatures. 
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Chapter 1 


Introduction 


Sulphur dioxide is released into the atmosphere through various sources such as coal- 
fired power plants, sulphuric acid plants, petroleum refineries, metal ore smelters, etc., 
the major source being coal-fired power plants. Among the various methods for control 
of sulphur oxides emissions cleaning up of flue gases received maximum attention. 

Prior to 1980, the removal of sulphur dioxide by absorption was usually carried out 
by wet scrubbing. Owing to a number of problems and disadvantages associated with 
wet scrubbing, attention has now shifted to dry scrubbing which minimises or eliminates 
scaling, corrosion and sludge handling problems a.sso dated with wet scrubbing and also 
eliminates both the reheating requirements and high-pressure drop conditions. 

The increased interest in dry flue gas desulphurisation (FGD) processes has led 
to the need for evaluating various sorbents. Evaluation of various sorbents under 
various process conditions may be facilitated if their performance can be predicted 
theoretically, using an appropriate model. Things are much more simplified if aid is 
taken of computer through the devlopment of a software package. 

The design of packed bed reactor involving gas-solid reactions poses considerable 
challenge because of the difficulties in characterising the transport and reacton pro- 
cesses of such reactors in a realistic manner . Szekely et al.(1976) quoted that in many 
instances even the individual steps involved in the overall desulphurisation process in 
a packed bed reactor are quite complex and inadequately understood. So the task of 
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presenting a workable model seems to be a formidable undertaking . Indeed, the state 
of affairs precludes the development of an entirely general ‘all purpose’ model which 
incorporates all the phenomena involved. 

In this work, a computer software has been developed to solve the equations repre- 
senting the grain model given by Szekely et al.(1976) and that has been used to predict 
the performance of the two sorbents ,namely, active sodium carbonate and calcined 
magnesium carbonate under various process conditions. A comparison has also been 
made between the predicted and experimentally observed values. 



Chapter 2 


Background Information 


2.1 Introduction 


Air pollution is woven throughout the fabric of our modern life. The major cause of 
all air pollution is combustion. Sulphur dioxide is one of the air pollutants produced 
due to the impurities present in the fuel. Recognition of sulphur dioxide as one of the 
pollutants is not a recent phenomenon. Episodes of severe sulphur dioxide pollution 
date back to several decades which led to efforts to control its emission. 


2.2 Harmful Effects of Sulphur Dioxide 


Various animal species, including man, respond to sulphur dioxide by bronchoconstric- 
tion, which may be assessed in terms of a slight increase in airway resistance. Sulphuric 
acid produced due to the presence of humidity is a much more potent irritant to man 
than sulphur dioxide. 

Different species and varieties of plants vary in their sensitivity to sulphur dioxide 
exposures. Alfalfa is the most sensitive. Entry through stomata causes injury to leaf 
and may lead to plant damage and loss of crop yield (Wark k, Warner, 1981). 
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Sulphur compounds are responsible for major damage to materials. Sulphur oxides 
accelerate metal corrosion, cause the loss of tensile strength of nylon and cotton fab- 
rics. Sulphurous or sulphuric acids are capable of attacking a wide variety of building 
materials including limestone, marble, roofing slabs and mortar. (Offen et. ah, 1987). 

Droplets of sulphuric acid mist formed by photochemical reactions between sulphur 
dioxide and other pollutants, contribute significantly to the reduction in visibility. 

Sulphuric acid formed in the atmosphere return to the ground as sulphates by means 
of dry and wet deposition. Acid deposition has a considerable impact not only on 
surface waters, but also on soils and ground water, mobilising micro nutrients and 
heavy metals and thereby impoverishing the status of the soils and poisoning drinking 
water. 


2.3 Sources of sulphur dioxide 


The sources of the oxides of sulphur are natural as well as anthropogenic, the latter 
being of prime importance to a pollution control engineer. The natural sources con- 
tributing to sulphur dioxide pollution are forest fires, volcanic eruptions and conversion 
of hydrogen sulphide into sulphur dioxide in the atmosphere. The past records and 
data indicate that nearly three fourth of anthropogenic sulphur dioxide emissions orig- 
inate from the combustion of fossil fuels (Wark and Warner, 1981). Sulphur is present 
in all natural mineral oils and coal with a composition varying from 0.1 to over 5 
percent. Fuel combustion in stationary sources and industrial processes are principal 
contributors of sulphur dioxide from human activities. About 69% of man-made sul- 
phur dioxide emissions in India are attributed to the thermal power generation during 
1990-91. Transportation contributes little to the antrhropogenic sulphur oxides in the 
atmosphere, because the sulphur in the gasoline is low. 


2.4 Flue gas desulphurisation systems 


Flue gas desulphurisation (FGD) is presently the most commonly used technology to 
comply with sulphur dioxide emission requirements. FGD processes may be grouped 
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according to two classifications (i) throwaway and regenerative and (ii) wet, wet/dry 
and dry. 

In throwaway processes, as the waste product formed is discarded, fresh chemicals 
must be continually added. In regenerative processes, removal agent can be continually 
regenerated in a closed-loop fashion. Throwaway processes have been found to be 
cheaper and less complex but involve wa.ste disposal problems. Regenerative processes 
have been commercially tested but have been found to be much costlier . Wet or dry 
processes are differentiated simply on the basis that whether the product formed is dry 
or wet. 


2.4.1 Wet scrubbing 


In this, removal of sulphur dioxide is achieved through chemical absorption by an al- 
kaline solution. Even though wet scrubbing technology has been commercially proven, 
corrosion, erosion, scaling, large water requirements and relatively higher costs are 
some of the problems and disadvantages associated with it. 


2.4.2 Wet/Dry scrubbing 


In this technology water requirement is reduced to one half of that for wet scrubbing 
with a view to reduce the problems associated with the use of circulating alkaline water 
solutions. In this, water in the solution sprayed into the flue gas stream evaporates by 
the heat of flue gases resulting in dry waste. 


2.4.3 Dry scrubbing 

Dry scrubbing processes carry spray drying philosophy one step further by eliminating 
water usage entirely. In dry scrubbing processes, a dry alkaline powder is injected 
directly into the flue gas stream. The alkaline particles react with sulphur dioxide 
while suspended in gas stream. Dry waste product from this reaction is , subsequently, 
collected in a particulate collection device and the clean flue gas is vented to the 
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atmosphere. 


2.5 Advantages of dry control technologies 


Although the reactivity of dry scrubbing agents will naturally be somewhat lower than 
the reactivity of wet sorbent solutions used in wet scrubbing, dry scrubbing olfers 
following advantages : 

(i) The absence of water results in the minimization or elimination of corrosion, 
erosion, plugging and scaling problems. 

(ii) Dry scrubbing processes control sulphur dioxide and particulates in a single piece 
of equipment. The capital costs of the process are projected to be 30 to 50% less than 
those of wet scrubber or spray drying systems. 

(iii) Dry scrubbing process eliminates both the reheating requirements and high 
pressure drop conditions. 


2.6 Dry desulphurisation in a packed bed 


Dry sulphur dioxide removal in a packed bed of solid reactant involves gas-solid re- 
action. Gaseous pollutant reacts with the solid reactant and ,thus, gets removed. A 
gas-solid reaction system may involve several intermediate steps. Typically, these in- 
termediate steps involve the following : 

REACTION- DIFFUSION 

(1) Gas phase mass transfer of the gaseous reactant from the bulk of the gas stream 
to the external surface of the solid particle. 

(2) (i) Diffusion of the gaseous reactant through the pores of the solid matrix, which 
could consist of a mixture of solid reactants and products. 
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(ii) Adsorption of the gaseous reactant on the surface of the solid matrix. 

(iii) Chemical reaction at the surface of the solid matrix. 

(iv) Desorption of the gaseous product from the surface of the solid matrix. 

(v) Diffusion of gaseous reaction product through the pores of the solid matrix. 

(3) Gas phase mass transfer of the gaseous product from the external surface of the 
solid to the bulk of the gas stream. 

HEAT TRANSFER 

(1) Convective heat transfer between the gas stream and the surface of the solid 
particle, and 

(2) Conduction heat transfer within the solid reactant product matrix. 
STRUCTURAL CHANGES 

The reaction and heat transfer processes could lead to structural changes, such as 
sintering or changes in the pore structure, which in turn could have a marked effect on 
the overall reaction rate. 

The mass transfer step has been extensively studied and is perhaps the best un- 
derstood. Pore diffusion is inherently much more complex than diffusion in liquids or 
gases and, as a consequence, is much less understood. Some of the complicating factors 
are immediately apparent : 

1) The volume occupied by the solids is not available for diffusive transfer. 

2) The actual diffusion path will not follow a straight line but will be quite tortuous, 
and the extent of this tortuosity will, necessarily, depend on the pore structure of the 
solids. 

The structural changes that may accompany gas-solid reactions can be quite complex 
and may be in the form of sintering, swelling, softening and cracking. Sintering implies 
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the phenomenon by which a porous compact increases its density as a consequence of 
being held at an elevated temperature, below its melting point. While a great deal of 
useful work is being done on sintering, the development of the field has not yet reached 
the stage where these models can be applied directly to the modelling of gas-solid 
reactions in which structural changes are important(Szekely et.al., 1976). 

Large number of experimental investigations on the process development with vari- 
ous types of reactive solids have been reported . However, very few studies deal with 
modelling and simulation of performance of a solid sorbent for desulphurisation in a 
packed bed. 

One of the common materials used for sulphur dioxide removal is lime in its various 
forms. Marrier and Dibles (1974) have reported some data on sulphation of calcium 
oxide and magnesium oxide. Kinetics of reaction of sulphur dioxide with calcined 
limestone has been studied by Borgwardt (1970) and that with calcium oxide has been 
studied by Borgwardt and Harvey (1972). 

A structural model was developed by Hartman and Coughlin (1974) to correlate 
the data on the reaction of sulphur dioxide with porous limestone in excess air at high 
temperature. Pigford and Sliger (1973) analysed the data of Borgwardt (1970) on the 
absorption of sulphur dioxide with porous calcium oxide and interpreted the data in 
terms of two diffusion processes - diffusion of sulphur dioxide through the pores and 
that through a layer of solid reaction product, coupled with chemical reaction. 

Hartman and Coughlin (1976) proposed a grain model incorporating the chemical 
reaction and diffusional resistances. Alvfors and Svedberg (1988) have further included 
the presence of inert matter in the particles. Ghosh et al. (1989) did modelling 
and simulation for desulphurisation of flue gas using dolomite. Dam Johansen et al. 
(1991) presented a modified grain model for the reaction between limestone, sulphur 
dioxide and oxygen taking into account the diffusion and external mass transfer to the 
particle surface. Karlegard and Bjerle (1994) compared the experimental results for 
desulphurisation of synthesis gas using zinc ferrite in apacked bed reactor, with those 
simulated on the basis of unreacted core model. 



Chapter 3 


Scope of the work 


Significance of modelling and simulation of performance of a solid sorbent for desul- 
phurisation needs no emphasis in view of the need for evaluating various sorbents in 
the dry flue gas desulphurisation systems. It ,also, helps in developing a better under- 
standing of the processes involved. 

Present work aims at predicting breakthrough curves for sulphur dioxide in a packed 
bed reactor of solid sorbents with the following objectives : 

(i) Developing a customized and interactive software package to solve the involved 
differential equation and obtain the breakthrough curves, 

(ii) To obtain predicted breakthrough curves for sulphur dioxide using the two solid 
sorbents, namely, active sodium carbonate and calcined magnesium carbonate under 
varying conditions of temperature and inlet sulphur dioxide concentration using the 
software developed and taking the basic parameters pertaining to the system like, rate 
constant, effective diffusivity etc. required for prediction from the literature (Bhaskar, 
1994) and 

(iii) Comparison of predicted breakthrough curves with the experimental ones, avail- 
able in the literature cited above. 
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Chapter 4 


Sulphation of solid sorbents and 
the packed bed model 


4.1 Sulphation of active sodium carbonate 


When sodium bicarbonate is exposed to flue gases at higher temperatures, thermal 
decomposition of the solid occurs producing the micrograin structured active sodium 
carbonate. Thermal decomposition of sodium bicarbonate above 115 °C is reasonably 
fast at atmospheric pressure yielding chemically active anhydrous sodium carbonate 
according to the equation (Bhaskar, 1994) 


2NaHC0z (s) — ^ NazCOz (s) + CO 2 (g) + H 2 O 


(4.1) 


AG = 6.79 kcal. 

The sodium carbonate formed according to above equation is known to have high 
surface area. 

Active sodium carbonate reacts readily with sulphur dioxide and oxygen yielding 
sodium sulphate according to the reaction 


10 



CHA PTER 4. SULPHATION OF SOLID SORBENTS AND THE PACKED BED MODELl 


Na^CO^ (s) + SO 2 (g) + ^-O, { 9 ) Na,SO^ (s) + CO 2 { 9 ) 


(4.2) 


AG = -76.45 kcal. 


4.2 Sulphation of calcined magnesium carbonate 


At higher temperatures magnesium carbonate decomposes according to the reaction 


MgCO^ (s) MgO (s) + CO 2 {g) (4.3) 

AG = 11.55 ibca/. 

In the presence of oxygen, the reaction of magnesium oxide with sulphur dioxide is 
given as 

MgO (s) + SOi (g) + jO, (<,) — » MgSO, (s) (4.4) 

AG= -71.95 kcal. 


4.3 Packed bed model 


Using a shell balance (Szekely et al., 1976), variation of concentration of the pollutant, 
A, in an isothermal packed bed with plug flow has. been shown as, 


dg ^ ^ ^ ^ ^ b at 


(4.5) 


where, 
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b stoichiometric coefficient 

Ca concentration of solute A in the fluid phase 

S cross section of the bed 

t residence time 

(7 superficial velocity of fluid through the bed 

X fractional conversion of solid 

y distance along the length of the packed bed 

Pb molar density of reactant solid 

Cp porosity of the pellet 

e„ porosity of the bed 

oj fractional volume of reactive solids in the bed 


In dimensionless form the above equation reduces to (Ghosh et ah, 1989) 


, a I A 

_ + a— + /3V. = 0 


(4.6) 


with the initial and boundary conditions 


0(0,0 = i = /(^*), 
'/’(?/*> 0) = 0 = 
where, 


( bkCA.AA 

\ ) 


(4.8) 


* 


A„ 


U 


iO 


F V 

A 3'' 9 / 


(4.9) 


y 


(1 -ep)(l -e)y 
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a is a distributed parameter and is equal to 


b{ep -f €.v)Ci 


w(l -Cp)(l - e„)pB 


(4.10) 


0 


ri 


^(1 - ( 2(1 - - 2 ) 

+(0.21 - 0.62X)(1 + a^)exp ( -0.9 (in^ 

\ \ 1.08 


-1 


(4.11) 


a = Ro 


(^ ^p) ^9 


11/2 


(4.12) 


The solution of Equation 4.6, through Laplace transformation has been obtained as 
follows, 


Cm dt‘ 


= 0 e 


-0yl 


(4.13) 


with the initial condition. 


X = 0 when t = 0 


(4.14) 


The above equation represents the packed bed model and its solution gives the variation 
of outlet concentration of gaseous reactant with time. 



Chapter 5 


Results and Discussion 


5.1 Salient features of the software package 


The software package developed is 
develop the breakthrough curves for 


an interactive, integrated and customized one to 
desulphurisation of flue gas in a packed bed reactor. 


Some of its salient features are : 


. It combines both the graphics as well as the problem solving features into one. 

. It prompts the user for the input at various points of its execution. 

. Input can be given either through the terminal or through input data-files. 

t • u fiTTip solution is 

• r ill - r ^ — 

as is required. 

. It has the option tor creating hardcopy of the output i.e the output can he printed 
on the plotter or the laser printer. 
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It creates *.plt file to be directly sent to the plotter for printing and *.bit file 
for the laser printer. The *.bit file has to be first translated using the pcltrans 
command before sending the request to the laser printer. 

Source code of the software developed is given in appendix A. 


5.2 Sample simulations 


Equations 4.13, 4.11, 4.12, 4.9 and 4.8 have been used for prediction of breaJk- 
through curves. These are given here, again, for their easy reference for discus- 
sion. 

Various basic parameters such as k,rate constant, Z>e, effective diffusivity, etc. 
appearing in these equations are given in appendix C. 


V’ 


where. 


tp is the dimensionless outlet concentration. 




^(1 - -f < 7 ^ ( 2(1 - - 2 ) 

.O 

+(0.21 - 0.62.Y)(1 + (7“)eip (-0.9 ) 


(5.1) 


(5.2) 


where, ,5 is a function of X, fractional conversion and a, dimensionless 

reaction modulous. 


cr 


Ro 


(i-g ]» 

\yj D 


(5.3) 


Vl = 



(1 ~ ^p) (1 - ^v) ^ 


(5.4) 


where. 


t/2 is the dimensionless lingth of the reactor. 
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5.2.1 Simulations for different values of a [y^ = 1) 

Predicted breakthrough curves for different values of cr are given in Table B.l 
and shown in Figure 5.1. It can be observed that breakthrough is delayed for 
lower values of cr and is quick for higher values. 

From the expression for cr and , it is clear that keeping constant, a can be 
made to vary with Ro, radius of the particle and De, effective diffusivity. Effects 
of Ro and Dg on breakthrough curves are discussed below. Explanations have 
been given both on the basis of concept as well as on the basis of mathematical 
equations representing the model. 

Effect of Ro 

Conceptually, for a given amount of solid reactant surface area available for re- 
action depends on the size of the particle. Surface area available being greater 
in the case of smaller particles and smaller in case of bigger particles. Hence, 
breakthrough will be delayed in the case of small particles and quick in the case 
of bigger particles. 

From the equations given in the proceeding section, it is clear that for large Ro , 
cr is large and hence /3 is small. Ca^ / would be greater for small value of 0. 
Hence, breakthrough would be quick in case of big particles and vice-versa. 

Effect of Dg 

Greater the effective diffusivity, more would be the access for the gas to the solid 
reactant and ,hence, more delayed would be the breakthrough curve. Reverse is 
the case for small effective diffusivity. 

Greater the value of De, smaller the value of a and ,hence, smaller the outlet 
concentration which means a delayed breakthrough curve. 

5.2.2 Simulations for different values oi y*^ { a = 1 ) 

Some predicted breakthrough curves for different values of are given in Table 
B.2 and shown in Figure 5.2. It is observed that breakthrough is more and more 
delayed with increasing t/^. 

From the expressions for a and y'l, it is clear that while keeping cr constant, f/£ 
can be made to vary with L j U , residence time, u;, fractional volume of reactant 
solid, and void fraction. 



Dimenaionleas Outlet Cone., C^£7< 
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IHrnerisumlesa Diatance, = 1.Q 



Sigma^O.Z 

Sifpna^O.4 

Sigma^0,8 

Sigmas f .2 

Sigmas f.$ 

SigmorS^Z^O 

— — Sigmaa^.O 

— *— Sigm4x:=^6.8 

— — Signw^fO.O 


Dimensionless Time, 


Figure 5.1: Predicted breakthrough curves for difFereut values of dimensionless reaction 
modulous, or 
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JHrneTiaionleas Reaction Modulus^ Sigma, = UO 
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Y£^ = 0,2 

n* = 0.4 

- - = 0.8 
yi* =s UZ 

n* = 1.6 

n* = 2.0 


Dimensionless Time, t* 


Figure 5.2: Predicted breakthrough curves for different values of dimensionless dis- 
tance, j/2 

Effect of residence time 

Greater the residence time, more would be the opportunity for the gas to re- 
main in contact and react with the reactant solid and hence lesser would be the 
outlet concentration. Therefore, breakthrough curve would be delayed for large 
residence time and vice-versa. 

From Equation 5.4, it is clear that greater the residence time, greater would be 
the value of From eqn. 5.1, outlet concentration would be smaller for greater 
yl and , hence, a delayed breakthrough curve. 

Effect of uj 

More the amount of reactant solid in the packed bed, more would be the oppor- 
tunity for the gas to react with it and result would be a delayed breakthrough 
curve. 

Greater the value of a;, greater would be the value of yl and , hence, a delayed 
breakthrough curve. 


Effect of 6,, 

If the void fraction is more, gas solid contact would be less and it will lead to 
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higher outlet concentration or a quick breakthrough. 

Greater the value of e„, lesser would be the value of ?/£ and hence a quick break- 
through. 


5.3 Comparison of predicted and experimental 
breakthrough curves for active sodium carbon- 
ate. 

5.3.1 At different temperatures 

Breakthrough curves are predicted at different temperatures 150, 200, 250, 300 
and 350 °C. Inlet SO 2 concentration and space velocity being 3500 ppm and 
900 hr~^, respectively. Predicted breakthrough curves at different temperatures 
are shown in the Figure 5.3. These have been compared with the experimentally 
observed ones. Comparison is given in the Tables D.1-D.5 and shown in the 
Figures 5. 4-5.8. 

Rate constant and effective diffusivity are the two parameters which vary with 
temperature. Variation of these parameters with temperature are given in Ap- 
pendix C. These values of rate constant and effective diffusivity were used for 
evaluating a and j/£ which in turn were used for predicting breakthrough curves. 
Effect of both the parameters is in the some direction. 

As expected, predicted breakthrough curves are delayed and less sharp at higher 
temperature (Figure 5.3). It can also be seen from the Figures 5. 4-5. 8 that the 
predicted breakthrough curves match fairly well with the experimentally obtained 
ones except for T = 350 °C, Departure of predicted breakthrough curves from 
the experimental one at 350 °C is attributed to the phenomenon of sintering. 
Sintering of pellets has been confirmed through scanning electron micrographs 
taken by Bhaskar (1994). Sintering results in a less porous and dense and compact 
solid reactant leading to poor diffusion into the reactant and ,hence, poor reaction 
of gas with it. 


5.3.2 At different Concentrations 

Breakthrough curves are predicted at different inlet SO 2 concentrations of 1500, 
3500, 6000, 1800 ppm and at a temperature and space velocity of 200 '^C and 
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Figure 5.3: Predicted breakthrough curves for sulphatiou of active sodium carbonate 
at different temperatures, (Ca, — 3500 ppm) 



Figure 5.4: Breakthrough curve for sulphation of sodium carbonate at 150 “’C. 
f Source : Bhaskar,1994. 
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Figure 5-5: Breakthrough curve for sulphation of sodium carbonate at 200 °C. 
f Source : Bhaskar,1994. 



cj 



Dimensionless Time, t* 


Figure 5.6: Breakthrough curve for sulphation of sodium carbonate at 250 '’C. 
f Source : Bhaskar,1994. 
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Figure 5.7: Breakthrough curve for sulphation of sodium carbonate at 300 °C. 
f Source : Bhaskar,1994. 



Dimensionless Time, t* 


Figure 5.8: Breakthrough curve for sulphation of sodium carbonate at 350 ‘‘C. 
f Source : Bha.skar,1994. 
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Figure 5.9: Breakthrough curves for sulphation of active sodium carbonate at different 
concentrations, {T = 200 ^C) 


900 respectively. Comparison of predicted and experimental breakthrough 
curves has been presented in Table D.6 and Figure 5.9. As expected, at higher 
values of inlet sulphur dioxide concentration, breakthrough has occurred faster. 
It is also seen that predicted curves match fairly well with the experimental ones. 


5.4 Comparison of predicted and experimen- 
tal breakthrough curves for calcined magnesium 
carbonate. 

5.4.1 At different temperatures 

Breakthrough curves have been predicted at different temperatures of 650, 750, 
and 850 °C. Concentration and space velocity are 9000 ppm and 600 Iir~^ , 
respectively. Predicted breakthrough curves at different temperatures are shown 
in Figure 5.10. 

Predicted and experimental breakthrough curves have been compared and pre- 
sented in Tables D.7-D.9 and Figures 5.11-5.13. 

As expected, removal of sulphur dioxide is better at higher temperatures owing 
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Figure 5.10: Breakthrough curves for sulphation of calcined magnesium carbonate at 
different temperatures, (Ca; = 9000 ppm) 

to higher reaction rates and higher effective diffusivities. Experimental and pre- 
dicted curves match fairly well except in the case of those at 850 °C. Again, the 
reason for this departure is sintering which has not been accounted for in the 
model owing to less understanding of the phenomenon. 

5.4.2 At different Concentrations 

Experimental and predicted breakthrough curves for sulphation of calcined mag- 
nesium carbonate at different concentrations of 2000, 9000, and 2200 ppm are 
presented in Table D.IO and Figure 5.14. Temperature and space velocity are 750 
°C and 600 hr~^, respectively. Predicted curves match fairly well with the ex- 
perimental ones and breakthrough is found to occur faster and sharper at higher 
inlet sulphur dioxide concentration. 
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Figure 5.11: Breakthrough curves for sulphation of calcined magnesium carbonate at 
650 °C. 

f Source : Bhaskar,1994. 



Figure 5.12: Breakthrough curves for sulphation of calcined magnesium carbonate at 
750 “C. 

t Source : Bhaskar,1994. 
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Dimensionless Time, i* 


Figure 5.13: Breakthrough curves for sulphation of calcined magnesium carbonate at 
850 “C. 

t Source : Bhaskar,1994. 
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Figure 5.14: Breakthrough curves for sulphation of calcined magnesium carbonate at 
different concentrations, {T = 750 ^ C) 

t Source : Bhaskar,1994, 
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5.5 Criticism of the present study 

As is clear from the Equation 5.1, dimensionless outlet concentration varies 
exponentially with ^ and Theoretically, at time t=0 dimensionless outlet 
concentration would be zero only if the product /JyJ, is infinity which is im- 
practical. So, dimensionless outlet concentration would only be nearly equal 
to zero at time t=0 as ^ or y£ or both become large which can be observed 
in the case of active sodium carbonate and sulphur dioxide system (espe- 
cially at higher tempyeratures) and not so in the case of calcined magnesium 
carbonate and sulphur dioxide system. 

From Equation 5.2 it is clear that ^ is inversely proportional to cr. As a 
increases, ^ would decrease and vice-versa. Therefore, breakthrough would 
be delayed and less steep for small a and vice-versa which is also supported 
by Figure 5.1. Then from Equation 5.3 , one would be tempted to say that 
since a is directly proportional to k, breakthrough would be delayed and less 
steep for smaller value of k. As is clear from the Equations 5.3 and 5.4, 
variation in k affects both a as well as y^. Effect of k on y£ is more than 
that on a because y2 varies with k raised to the power 1 and cr varies with 
k raised to the power |. Therefore, effect of y^ due to variation of k on 
breakthrough curve would dominate that due to cr and, hence, breakthrough 
would be delayed and less steep for higher values of k and vice-versa. For 
other parameters (not expained in the sections 5.2 and 5.3) the breakthrough 
curves’ trend can be explained on similar lines. 

Figure 5.2 shows that breakthrough is quicker for smaller value of yl or 
smaller value of k which is again supported by both the Equations 5.1 and 
5.4. 

Departure of simulated and experimental breakthrough curves may also 
be attributed to, besides sintering, change in porosity due to difference in 
molar volumes of the reactant and product (especially in the case of cal- 
cined magnesium carbonate and sulphur dioxide system), pore closure and 
changing effective diffusivity. 



Chapter 6 


Conclusions and Suggestions for 
Future Work 


On the basis of the discussion presented in the last chapter it may be concluded 
that the model presented here for simulation of breakthrough curves in a packed 
bed reactor simulated fairly well the experimental results for the two solid reac- 
tants namely, active sodium carbonate and calcined magnesium carbonate. The 
observed departure of simulated results from the experimental data , especially, 
at higher temperatures in the case of both the sorbents may be attributed to the 
structural changes such as sintering, change in porosity etc. occuring in the solid 
reactant during the course of the reaction. 


A more comprehensive model which takes into account the non-isothermal con- 
ditions and structural changes should be developed for better simulation. 


With the help of computer, a number of sets of curves may be obtained by varying 
different parameters involved one by one. Then to get predicted breakthrough 
curves for a specific reactant and under different process conditions, interpolation 
may be resorted to. 
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Appendix A 

Source Code of the Software 
Package 


tinclude <starbase.c.h> 

#iiiclude <math.h> 

#include <stdio.h> 

#include <stdlib.li> 

#iiiclude <string.h.> 

#defiiie Edge TRUE 

tdefine SimplexSansSerif 4 

#defiiie BoldSaasSerif 6 

#defi]ie EndOfLine FALSE 


struct s 

1 { char 

♦fntl , 

♦fnt2, 

♦fnd; 

} filename [10] = { 

■["tl.l". 

"tl.2". 

"dl.o"}. 

■C"t2.1", 

"t2.2". 

"d2.o">. 

{"t3.1", 

"t3.2". 

"d3.o"}, 

{"t4.1". 

"t4.2". 

"d4.o"}. 

{"t5.1", 

"t5.2". 

"d5.o"}. 

{"te.l". 

"t6.2", 

"d6.o"}, 

■C"t7.1", 

"t7.2", 

"d7.o"}. 

{"tS.l", 

"t8.2". 

"dS.o"}, 

{"t9.1". 

"t9.2". 

"d9.o"}. 


{"tlO. 

1", "tlO.2", 


"dlO.o"} 

}; 


FILE *ptdatain,*ptdatainl,*ptdataout [20] ,*pttextoTitl [20] , 
*pttextout2[20] , *original, ♦datafile [10] ; 
char xstringl[50]="DIMENSI0NLESS TIME t+". 
xstring2[50]="TIME (min.)", 

ystringl [50]="DIMENSI0NLESS OUTLET CONG. Co/Ci", 
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ystring2[50]="FRACTI0NAL CONVERSION X", 
mainstringl [50] =" BREAKTHROUGH CURVES" ; 
mainO 
{ 

12345678901234567890123456789012345678901234567890123456789012345 

/* priiitf("\ii\ii\t\tDEVELOPED BY\n\n\ii") ; 

systemCbanner SATISH KUMAR") ; */ 

system("clear") ; 

prt() ; 

priiitf("ENTER\ii\t\'l\’\t FOR SOLVING DIFF. EQN. ALS0\n\t\'2\'\t 
FOR DRAWING GRAPHS 0NLY\n ? "); 
scanf ("*/,c" j&rrnod) ; 
scanf ("y.c" , Mummy) ; 
if Crmod=='l’) 

{ 

priixtf("\iiENTER\ii\t\'l\'\tIF WANT TO HAVE ONLY THE SOLUTION OF 
DIFF. EQN. F0R\ii\ \rTHE GIVEN SETS OF DATA( ENTER 1 IF 
NO. OF SETS OF DATA IS MORE THAN FIVE )\n\n\ \r\t\'2\'\t 
IF WANT TO SEE GRAPH ALSO\n\n ? "); 
scanf ("%c" ,4r) ; 
scanf ("Xc" .Mummy) ; 

/♦ printf("\n\nENTER\n\t\'l'\tFOR EXECUTING VERSION OF THE MODEL 
IN WHICH SIGMA IS KEPT CONSTANT\n\ \t\’2\^\tF0R EXECUT 
ING THE VERSION IN WHICH SIGMA VARIES WITH X.\n ? ") ; 
scanf (" y.c" ,&rmodel) ; 
scanf ("5ic" .Mummy) ; 
if (irmodel== ’!’) */ 

difflO; 

/* else diff2(); */ 
if (r=='lO 
exit(-l) ; 

} 

/*♦♦♦**♦♦*♦♦*♦*♦ Input of main.x-axis.y-axis titles 

printf ("\nHERE IS THE LIST OF FILES IN THE PRESENT DIRECTORY 
\n\n"); 
system("ls") ; 

printf("\n\nENTER THE NO. OF DATA-FILES ( 0,1,2,... )\n 
(excluding those generated by the slon. of diff. eqn.) 

: "); 

scanf ("y.d",&nf); 
scanf ('"/.c" .Mummy) ; 
if (rmod=®=’ lO 
tn = nf + ns; 

else tn=nf; 
for (j=n3; j<tn; j++) 

printf ("\nDATAFILE[y.d] = ", j+1); 
gets (f name [j] ) ; 

> 
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printf("\nENTER THE MAIN TITLE : "); 

printf("\nENTER\n\t\'l\'\tFOR \ "BREAKTHROUGH CURVES\"\n\t\' 
2\'\tF0R ENTRING SOME OTHER\n ? "); 
scanf ("’/.c" ,&r) ; 
scanf ('"/.c" ,&diinmiy) ; 
if(r==’lO 

strcpy(maiii_title,mainstringl) ; 
if (r=='20 

gets(main_title) ; 

printfC'XnENTER THE X-AXIS TITLE : ") ; 

printf ("\nENTER\ii\t\'l\'\tFOR \" DIMENSIONLESS TIME t* \"\ii 
\t\'2\’\tF0R \"TIME (min) ->\"\ii\t\’3\'\tF0R ENTRING 
SOME OTHERXn ? "); 
scanf ("7, c" ,&r) ; 
scanf ("7c" , Mummy) ; 
if (r=='lO 

strcpy(x_title,xstringl) ; 
else if(r=='20 

strcpy(x_title,xstring2) ; 
else 

gets(x_title) ; 

printf ("\nENTER THE Y-AXIS TITLE : ") ; 

printf ("\nENTER\n\t\'l\'\tFOR \"DIMENSIONLESS OUTLET CONG. 

Co/Ci \"\n\t\'2\’\tF0R \"FRACTIONAL CONVERSION X\"\n\t\ 
'3\'\tF0R ENTRING SOME 0THER\n ? ") ; 
scauif ("*/,c" ,&r) ; 
scanf ("7c" , Mummy) ; 
if (r=='lO 

strcpy(y_title,ystringl) ; 
else if(r=='2') 

strcpy(y_title,ystring2) ; 
else 

gets(y_title) ; 
f or(j=0; j<tn,* j++) 

{ 

printf ("\nENTER THE TITLE FOR THE 7dst CURVE\n (DEFAULT IS THE 
FILENAME) : ",j+l); 
gets (legend Cj]) ; 
if (strlen(legendCj] ) <= 2) 

/* legendCj] = fnameCj] ; */ 
strcpy(legendCj] ,fnameCj] ) ; 

} 

printf ("\nDO YOU WANT TO ENTER THE VALUE OF THE DIVIDER?\n\ 
XrEFAULT VALUE OF DIVIDER IS 100. RUN THE PROGRAMXnX 
XrONCE WITH DEFAULT VALUE. IF THE SIZE OF TICS ISXnX 
XrLARGE,NEXT TIME, ENTER A VALUE > 100 FOR THE DIVIDERXnX 
Xr AND VICE VERSA.XnXn(y/n) : "); 
scanf ("7c" ,&r_divider) ; 
scanf ("7c" , Mummy) ; 
if (r_divider == 'y') 

■C printf ("XnENTER THE VALUE OF DIVEDER : "); 
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scaaifC" •/.(!", Mivider) ; 
scanf j&dummy) ; 


/♦♦***♦♦♦*♦♦***♦♦♦ setting parameters **************************/ 

for (m=0 ;m<tn;m++) 

{ 

i=l; 

k=l; 

datafilelm] =f open (f name [m] ,"r") ; 
fscanf (datafileW ,"y.f’/,f",&X,&Y) ; 
x_data[m] C0]=X; 
y_dataCm] C0]=Y; 

while (f scant (datafile Cm] , '"/.f/.f", &X, &Y) != EOF) 

{ 

x_data[m] [i++]=X; 
y_dataCm] Ck++]=Y; 

} 

N Cm] =i ; 

xmax=x_dataCm] CO] ; 
ymax=y_dataCiii] CO] ; 
xmin=x_dataCni] CO] ; 
yiiiin=y_dataCm] CO] ; 
for(i=l,j=l; i<NCm] ; i++,j++) 

{ 

if (x_dataCm] Ci] > xmax) xm 2 ix=x_dataCm] Ci] ; 
if (y.dataCm] Cj] > ymax) ymax=y_dataCm] Cj] ; 

> 

local_x_max Cm] =xmax ; 
local_y_max Cm] =ymax ; 
for(i=l, j=l,* i<NCm] ; i++,j++) 

{ 

if (x.dataCm] Ci] < xmin) xmin=x_dataCm] Ci] ; 
if (y.dataCm] Cj] < ymin) ymin=y_dataCm] Cj] ; 

} 

local_y_min Cm] =ymin ; 
local_x_min Cm] =xmin ; 
fclose(dataf ileCm ] ) ; 

y 

f or(i=0 ; i<tn; i++) 

printf("\nNUMBER OF DATA IN THE FILE ’/.s = '/.d" ,fnameCi] ,NCi] ) ; 

xmax=local_x_max Co] ; 

ymax=local_y_maxCO] ; 

xmin=local_x_minCO] ; 

ymin=local_y_minCO] ; 

f or(i=l , j=l ; i<tn; i++, j++) 
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if (local_x_m 2 ix[i] > xmax) 
if (local_y_majcCj] > ymajc) 

abs_x_max = xmax; 
a-bs_y_max = ymax; 
for (i=l, j=l;i<tn;i++,j++) 

if (local _x_min[i3 < xmin) 
if (local_y_miii[j] < ymin) 

> 

abs_x_min = xmin; 
abs_y_miii = ymin; 


xmajc=local_x_maLx[i] ; 
ymax=local_y_maxCj] ; 


xmin=local_x_min[i] ; 
ymin=local_y_min[j] ; 




delta_x = (abs_x_max - abs_x_min) /divider; 
delta_y = (abs_y_max - abs_y_min) /divider; 
x_ticsize = delta.y; 
y_ticsize = delta_x; 

printfC"\n\nPRESENT CO-ORDINATE LIMITS ARE : \n\n"); 
printf ("\n\nABS-X-MIN = ‘/.fXnABS-Y-MIN = '/.f \nABS-X-MAX = ‘/.fXnABS-Y 
-MAX = ‘/f\ii\n\ ",abs_x_min,abs_y_min,abs_x_max,abs_y_m£ix) ; 
printf ("\nDO YOU WANT TO CHANGE THESE LIMITS ? (y/n) : "): 
scanf ("*/c" ,&r) ; 
scanf ("%c" .Adummy) ; 
if (r=='yO 

printf ("\n\nX-MIN = ") ; 
scanf ("i(f" ,&abs_x_min) ; 
scanf ("%c" , Mummy) ; 
printf ("\nY-MIN = ") ; 
scanf (“Xf" ,&abs_y_min) ; 
scanf ("%c" , Mummy) ; 
printf ("\nX-MAX = ") ; 
scanf ("'/f" ,&abs_x_max) ; 
scanf ("Xc" .Mximmy) ; 
printf ("\nY-MAX = "); 
scanf ("*/,f" ,&abs_y_maLx) ; 
scanf ("Xc" , Mummy) ; 

> 

if ((abs_x_max==abs_x_min) II (abs_y_max==abs_y_min)) 

{ printf ("\n\nmax value = min value\n\n") ; 
exit(-l) ; 

} 

/*for(j=0; j<=n; j++) 

printf C'xCXd] = Xf yCXd] * Xf \n" , j ,x_data[j] , j ,y_data[j] ) ; */ 
/♦*♦♦♦♦**♦+♦**♦♦***♦♦♦♦♦♦*****♦♦♦♦*♦*♦**♦♦♦**♦*♦*♦****♦**♦******/ 


printf ("\' ENTER THE DEVICE TYPE\' ?\n\nENTER\n\t\'ws\'\tFOR 

W0RKSTATI0N\n\ \r\t\'srx\'\tFOR \’SRX\’\n\t\'t3rx\'\tF0R 
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\^TSRX\'\]i\n\t? "); 
gets(r_device) ; 
if (strcmp(r_device, "ws")==0) 
r_d= ‘ 1 ’ ; 

else if (strcmp (r_device , " srx") ==0) 
r_d= ’ 2 ’ ; 

else if (strcmp(r_device, "tsrx")==0) 
r_d='3' ; 

else printf ("\n\nERROR IN INITIALISING DEVICE TYPEXnXn") ; 
switch.(r_d) 

{ 

case '1' : 
term_type = ws; 
printf ("TERM = y.s\n",ws); 
break; 
case '2' : 
term_type = srx; 
printf ("TERM = y.s\n",srx); 
break; 
case '3' : 
term_type = tsrx; 
printf ("TERM = y.s\n",tsrx); 
break; 

} 

Input of main, X- axis, y-axis titles ♦*♦****♦*♦♦/ 

printf ("\nWANT TO VIEW THE GRAPH\nOR\nWANT TO CREATE THE PLOT 
FILE\ TO BE PLOTTED ON THE PLOTTER \n\nENTER\n\n\t\ ' l\ ' 

\tFQR\ VIEWING THE GRAPH FIRST AND THEN\n\t\t CREATING THE 
PLOT FILE\ \n\n\t\'2\'\tF0R CREATING THE PLOT FILE ONLY. 
\n\n\t? "); 
rl=getchar() ; 
scanf ("’/.c" ,&dTimmy) ; 
printf ("\nWAIT !\n"); 
if(rl=='l’ && ag==’n') 

{ 

printf ("\n\' DO YOU WANT TO HAVE THE BITMAP FILE FOR PRINTING\' 
(y/n) : "); 
scanf ("Xc" ,&r_bit) ; 
scanf ("Xc" ,&dummy) ; 
if (r_bit=='y') 

printf("\n\n\’ ENTER THE NAME OF THE BITMAP-FILE\ ' : ") ; 
scanf ("Xs" , bit .filename) ; 
scanf ("'/.c" ,&dtimmy) ; 

lenb=strien (bit .filename) ; 
for(i=0;i<5;i++) 
bit.f ilenameClenb+i] = bitCi]; 
printf ("\nWAIT !\n"); 

} 
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else 

printfC"\nWAIT 

} 

again : if(ag=='y’) 

{ 

rl='2’; 

} 

if (rl=='2’) 

{ 

printf("\n\nENTER THE NAME OF THE PLOT FILE ; "); 

sccinf ("*/,s" ,plot_f ilename) ; 

lenp=strlen (plot .filename) ; 

scanf ("*/c" ,&dnmmy) ; 

f or(i=0 ; i<5 ; i++) 

plot _f ilename [lenp+i] = plt[i]; 

printf ( " \nWAIT ! \n" ) ; 

} 

switch(rl) 
case ‘ 1 ’ : 

printf ("\331i\33j ") , fflush.(stdout) ; 

if ( (f ildes=gopen ( "/ dev/crt " , OUTDEV , term.type , INIT) )== 
exit(-l) ; 
system("clear") ; 
break; 
case ' 2 ' : 
system("clear") ; 

printf ("\tWAIT !\n'’); 

if ( (f ildes=gopen (plot _f ilen 2 une , OUTDEV , "hpgl " , SPOOLED ) ) == 
exit(-l); 

break; 

} 

buffer_mode(f ildes,TRUE) ; 

/*♦*♦**♦*+♦**♦ SELECTING VIEWPORT, WINDOW ETC. ♦♦*♦***♦**++***/ 

vdc.ext ent (f ildes ,0.0,0.0,0.0,1.25,1.0,0.0); 
view.port (f ildes , 0 . 35* 1.25,0.35,0. 95* 1.25,0.85); 
view_window(f ildes, abs_x_min,abs_y_min,abs_x_max,abs_y_max) ; 
int erior.style (f ildes , INT.HOLLOW , Edge) ; 

rectangle(f ildes , abs_x_min , abs_y_min , abs_x_max , abs.y.maix) ; 

PLOTTING GRAPH *♦*♦***♦♦♦♦*♦♦***♦*♦♦***♦*♦****♦/ 

k=0; 

i=0; 

for (j=0;j<tn;j++) 

printf ("\n\nCOLOUR INDEX CAN HAVE VALUE RANGING FROM 0-256. \n\ 
COLOURS AVAILABLE ON THE PLOTTER ARE FROM 0-7.\n\n\ 
INDEX\n\ 


COLOUR 
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\n\ 

WHITE l\n\ 

RED 2\ii\ 

YELLOW 3\n\ 

GREEN 4\n\ 

CYAN 5\n\ 

BLUE 6\ii\ 

MAGENTA 7"); 

printfC'XnXnCQLOUR INDEX FOR THE FILE ’/.s : " ,fnameCj] ) ; 

scauf ("*/,d" ,&colour_ index) ; */ 

dataf ileCj] =f open(fnameCj] ,"r") ; 

fscanf (datafile [j] ,"'/,f'/.f " ,&X,&Y) ; 

move2d(f ildes,X,Y) ; 

whileCfscanf (datafileCj] , "*/,f'/,f", &X, &Y) != EOF) 

ifCj==6.) k=0; 
if(j==6.) i=0; 
lin.e_type(f ildes , i) ; 
line_color_iiidex(fiides,2+k) ; 
draw2d(f ildes, X,Y); 

} 

i++; 

k++; 

fclose(datafileCj]) ; 

> 

/*♦*♦♦*♦*♦*♦♦♦*♦♦ do the main, title ****************************/ 

interior_style(f ildes , INT.SOLID , Edge) ; 
text_font_index(f ildes , SimplexSamsSerif ) ; 
text_alignment(f ildes, TA_CENTER,TA_ CENTER.O . ,0 .) ; 
charact er_height (f ildes , . 03) ; 
text_color_index(f ildes,4) ; 

text2d(f ildes, 0.65*1. 25, .93,main_title,VDC_TEXT,End0fLine) ; 

/*♦*+**♦*♦♦***♦*♦♦ (Jo "the x~axis title **************>^:*********/ 

text _font_index(f ildes, SimplexSansSerif) ; 
text_alignment(f ildes, TA_CENTER, TA_CENTER,0 . ,0. ) ; 
charact er_height (f ildes , . 03) ; 
text_color_index(fildes,4) ,* 

t ext 2d (f ildes , 0 . 65* 1 . 25 , 0 . 3 , x_t itle , VDC.TEXT , EndOf Line) ; 

/♦♦***♦♦♦**♦*♦** (Jo the y-axis title ♦*♦♦♦♦*♦♦*♦**♦*♦♦*****♦***/ 

t ext _font_ index (f ildes, SimplexSansSerif) ; 
text_ 2 ilignment(f ildes, TA_CAP,TA_CENTER,0. ,0.) ; 
charact er.height (f ildes , . 03) ; 

text_orientation2d(f ildes, -1.0 ,0. ,0. ,1.0); 
text color_index(f ildes ,4) ; 

text 2 d(f ildes , 0 . 27*1 . 2 , 0 . 6 , y.title, VDC.TEXT , EndOf Line) ; 
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/:|c*3(c])c>f:*:(c*>tc******** TIC MARKS ********************************/ 

line_color_index(f ildes,4) ; 
clip.indicator (f ildes , CLIP.TO .VIEWPORT) ; 

Tics(f ildes, 'X’ ,abs_y_miix,abs_x_min,abs_x_mcLX,delta_x,x_ticsize, 10) ; 
Tics(fildes, 'X' ,abs_y_max,abs_x_min,abs_x_max,delta_x,x_ticsize, 10) ; 
Tics(fildes, ’Y‘ ,abs_x_min,abs_y_min,abs_y_m 2 LX,delta_y,y_ticsize, 10) ; 
Tics (f ildes ,’Y’ , abs.x.max , abs.y.min, abs.y.max , delt a.y , y.ticsize , 10) ; 

label the X tic labels ************************/ 

clip.indicator (f ildes, CLIP.TO.VDC); 
text .orientation2d(f ildes ,0.0, 1.0, 1.0, 0.0); 
chaxacter.height(f ildes, 4 * delta.y ); 
character.width(f ildes, 1.5 * delta.x) ; 
text_color.index(fildes,4) ; 

t ext _al ignment (f ildes , TA.CENTER , TA.TOP , 0 . 0 , 0 . 0 ) ; 
for(X=abs.x.min;X<=abs_x.max;X+=10.0 * delta.x) 

■C 

if (abs.x.max<=10 . ) 
sprintf (Strng, "'/,5.2f",X); 
else if (abs.x_majc<=40 . ) 
sprintf (Stmg, "*/,4. If" ,X) ; 
else 

sprintf (Strng, "%4.0f",X); 

text 2d (f ildes , X , abs.y.min , Strng , WORLD.COORDINATE.TEXT , EndOf 
Line) ; 

} 

sprintf (Strng, "7,5. If", abs.x.max); 

t ext 2d(f ildes, abs.x.mELX, abs.y.min, Stmg, WORLD.COORDINATE.TEXT 
, EndOf Line) ; 

/♦♦♦♦♦♦♦♦♦♦♦♦♦label the Y tic labels ***************************/ 

text.alignment(f ildes, TA.RIGHT,TA.HALF,0 .0,0.0) ; 
text. color. index (f ildes, 4) ; 

f or(Y=abs.y.min;Y<=abs.y_max + delta.y ; Y+=10 . 0 * delta.y) 

if (abs_y_max<=10 . ) 
sprintf (Strng, "'/.5.2f ",Y); 
else if (abs.y.max<=40.) 
sprintf (Stmg, "7.4. If " ,Y) ; 

oX S6 

sprintf (stmg, "%4.0f",Y); 

text2d(f ildes, (abs_x.min-2.*delta.x) ,Y,Stmg,WORLD_COORDINATE 
TEXT, EndOf Line) ; 

} 

flush.buffer(f ildes) ; 

/♦♦♦♦♦♦♦♦♦♦♦♦♦ viewing graph ♦♦★♦*******>i<***********************/ 
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/* printf ("PRESS ENTER TO VIEW THE GRAPH\n"); 
scanf ("y,c" , Mummy) ; 

printf ("\nD0 YOU WANT TO CLEAR THE SCREEN BEFORE VIEWING GRAPH 
(y/n) : 

scauif ('"/.c" ,&r) ; 
if (r==’y') 
systemC'clear") ; */ 

/****♦♦***♦♦ SECOND VIEW PORT FOR LEGEND **********************/ 

vdc_extent(fildes, .0, .0, .0,1.25,1. , .0) ; 
view.port (f ildes , 0 . 4*1 .25,0.05,0. 98*1 .25,0. 25) ; 
view_window(f ildes ,.0,. 0,100. ,100.); 
interior_style(f ildes, INT.HOLLOW, Edge) ; 
text_font_index(f ildes, SimplexSansSerif) ; 
text .alignment (f ildes, TA.LEFT, TA_CENTER,0 . ,0.) ; 
chaxacter.height (f ildes , 10 . ) ; 
cli.aracter_width(f ildes ,1.8); 
k=0; 
i=0; 

f or(j=0; j<tn; j++) 

i 

if(j>13) 
break; 
if (j ==6) 

k=0; 

i=0; 

} 

if (j>=6) 

{ 

line_type(f ildes, i) ; 
line_color_index(f ildes, k+2) ; 
move2d(f ildes,60 .0,90. “(5 .+12 .*i)) ; 
draw2d(fildes,70.0,90.-(5.+12.*i)) ; 
text_color_index(fildes,2) ; 

text2d(f ildes ,75 . 0 , 90 . - (5 . +12 . *i) , legend [j] , WORLD.COORDINATE 
.TEXT, EndOf Line) ; 
k++; 
i++; 

continue; 

} 

line.type(f ildes, i) ; 
line.color.indexCf ildes, k+2) ; 
move2d(f ildes , 5 . 0 , 90 . - (5 . +12 . *i) ) ; 
draw2d(f ildes , 20 . 0 , 90 . - (5 . +12 . *i) ) ; 
text_color_index(fildes,2) ; 

text2d(f ildes , 25 . 0 , 90 . - (5 . +12 . *i) , legend [j] , WORLD.COORDINATE 
.TEXT, EndOf Line) ; 
k++; 



APPENDIX A. SOURCE CODE OF THE SOFTWARE PACKAGE 


40 


i++; 

} 

make_picture_current(f ildes) ; 
f lush._buffer(f ildes) ; 

Th-ird view port ♦♦♦♦♦♦♦♦jk+^+^k*************/ 

vdc.extent (fildes, .0, .0, .0,1 .25,1. , .0) ; 
view.port (f ildes ,0.4*1.25,0.0,0. 98*1 . 25,0. 05) ; 
view_window(f ildes, .0, .0,60. ,5.) ; 
interior_style(f ildes , INT_H0LL0W,Edge) ; 
text_foiit_index(fildes,SimplexSaiisSerif) ; 
text_aligiment(f ildes, TA_ CENTER, TA_CENTER,0 . ,0 . ) ; 
ch.axacter_height(fildes,2.75) ; 
chaxacter_widtli(f ildes ,1.2); 
t ext _color_ index (f ildes, 2) ; 

text 2d (f ildes, 30. 0,2. 5, "DEVELOPED BY SATISH KUMAR" ,W0RLD_C0QRDIN 
ATE_TEXT,EiidOfLine) ; 

make_picture_current(f ildes) ; 

/:tc3tc4:******************« CREATING BITMAP FILE ******************/ 
if(rl=='l' ag=='ii' kk r_bit=='y') 

■C 

printf("\nWAIT !\n"); 
bitmap_to_file(f ildes, 

TRUE, 

0 , 

0 , 

bit_f ilename, 

TRUE, 

.25,1. ,820,1000, 

TRUE) ; 

} 

/*♦**♦*♦♦*♦♦♦♦**♦♦ CLEARING SCREEN FOR THE IMAGE ♦*************\ 

/* system("cle3Lr") ;*/ 

/* make_picture_curren.t (f ildes) ; */ 

/*sleep(2); */ 

/♦♦♦♦♦♦*♦♦♦♦* retrieving image from the file ♦*♦♦♦♦=«♦♦♦*♦♦♦**♦*/ 

/* file_to_bitmap(f ildes,*/ 

/* TRUE, */ 

/* 0 , 0 , */ 

/* bit.filenaune, */ 

/* 0 . 0 , 1 . 0 , */ 

/* TRUE); */ 

/♦♦♦»<*♦♦***** CLOSING THE DEVICE ♦♦♦*♦♦♦♦**♦♦*******♦**********/ 



APPENDIX A. SOURCE CODE OF THE SOFTWARE PACKAGE 


41 


if (rl==^lO 

if (gclose(f ildes)==0) 

priatfC'file closed successfully. Xu") ; 
else 

printf ("Error in closing file.Xn"); 
printfC'XnDO YOU WANT TO CREATE PLOT FILE ? (y/n) : ") ; 
ag=getchax() ; 
scanf ("‘/,c" j&dummy ) ; 
if(ag=='y') 
goto again; 

> 

if (ag=='y') 

if (gclose(f ildes)==0) 

printf ("file closed successfully .\n") ; 
else 

printf ("Error in closing file.Xn"); 
printf ("XnXnXnX' END OF THE PROGRAMMEX'XnXnXn") ; 

> 

Z********** SUBROUTINE Tics nc**^c****************************/ 

Tics (fildes, Axis .Location, Start , End, Belt a, Size, Major) 
int fildes; 
char Axis; 
float Location, 

Start , 

End, 

Delta, 

Size; 

int Major; 

#define Len (Tic'/,Major!=0 ? Size : 2.0*Size) 

#define move 0.0 

.#define draw 1.0 

{ 

int Tic; 
float X,Y; 
float TicMksCl200] ; 
float *TicPtr; 

TicPtr=*TicMks ; 
switch (Axis) { 
case'X’ ; case'x’ : 

f or (X=Start ,Tic=0 ; X<=End; X+=Delt a. Tic++) 

*TicPtr++ =X ,*TicPtr++=Location-Len . *TicPtr++=move ; 
*TicPtr++ =X , *TicPtr++=Location+Len , ♦TicPt r++=draw ; 

> 

break ; 

case 'Y' : case'y' : 

for(Y=Start,Tic=0;Y<=End;Y+=Delta,Tic++) 
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{ 

*TicP'tr++ =Locatioii-Len, *TicPtr++=Y, *TicPtr++=iaove ; 
*TicPtr++ =Locatioii+Len,*TicPtr++=Y,*TicPtr++=draw: 
} 

break; 

} 

polyliiie2d(fildes,TicMks,(TicPtr-TicMks) / 3, TRUE); 


void prt() 

{ 

char dummy; 
printf ('* 

\37***5»c:4c5|e3*<5»cs»c+%**a»c**5|c3|cjk****************a|c****************5|c********* 

\r************************************************************** 

\j‘*3(ea|e5fc5|c:4c a|cjf:*s»c5fc* 

\r****** name of the STUDENT : SATISH KUMAR 

\j^a*e:|e^:4e5»c3|e :te3(C3te)ie*:fe 

\r****** THESIS SUPERVISOR : Dr. D.K. GHOSH 

\l7**:*C3*e5tca»e 5|c:^c*5*c:*e:»e 

\2::fc3«c:|e:)|c:tc}fc;4c*3(e)|e3te*))c3(c:»c:4e3|c3)c)^3fe;4c3|e4::4e3|c:4(**3(e:fe:4e****3ie:>ie3ic*:ic3|e*9te:fc::{c:4e*^:4c**:>|c:4c**3tc3ie3fe^:4c** 
3(e 3fe :fc 3)c afc :(( :4e a)e 3(e ofe 3|c 9tc^ 4c 9|e )fe % :jc 9(c :{e :4c ^ 

♦ ♦♦♦*♦*♦* 4 c 4 c ♦♦♦* 4 e \ n \ 

4 c 4 e 4 c 4 c 4 c 4 c\n\ 

♦** TITLE ***\n\ 

4 c 4 t 4 c ***\ll\ 

\ j- 4c * 4c 4c 4c 4c 4c 4c 4c 4c 4e 4e 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4e 4c 4c 4e 4c 4c 4c 4c 4c 4c * 4c 4e 4c * 4c 4e 4c 4c 4e 4c ♦ 4c 4c 4c 4e 4c 4c 4c * 4c 4c 4c 4c * 4c 4c 4c 4e 4c 
\ j* 4e 4c 4c 4c 4c 4c 4c 4c * 4c 4e 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4e 4c 4c 4c 4c 4c 4c 4c 4c 4e 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 

4c 4c 

SIMULATION OF BREAKTHROUGH CURVES FOR DRY *♦ 

4c 4c 

DESULPHURIZATION OF FLUE GAS USING * 

4c 

CALCINED MAGNESIUM CARBONATE ** 

4c 4c 

\j^4e4e4c4c4e4c4c4e4c4c4c4c4c4c4c4c4c4c4c4c4e4c4c4c4c4c4c4c4e4c4c4c4c4e4c4c4c4c4c4c4c4c4c4e4c4c4c4e4e4c4c4c4:4c4c4e4c4c4e4e4c4c 
\l^4c 4c 4c 4c 4c 4c 4c 4c 4C 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4e 4c 4c 4c 4c 4t 4e 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c 4c 4c4c4e 4c 4c 4c 4c 4c 4c 4t 4c 4c 4c 4c 

printf ("XnPress ENTER ..An**); 
scanf (’*y,c“ , &duinmy) ; 
system('*cleax**) ; 
printf (’* 

} 

/ 4 c 4 c 4 c 4 c 4 c 4 c 4 c 4 c 4 c 4 c 4 c 4 c 4 c 4c 4 c 4e 4c 4c 4c 4c4c 4c 4c 4c 4c 4c 4c 4c 4c4e 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c 4c 4c 4c 4c 4c4c 4c 4c 4c 4c 4c 4c/ 


\r** 

\r** 

\r 4 c 4 c 

\r** 

\j-4c4c 

\r** 

\r4c4c 


void difflO 

{ 

float belO; 
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float fnl() ; 
void solvelO; 

printf ("\nENTER\n\t\'l\'\tIF YOU WANT TO ENTER DATA FROM THE 
TERMINAL\n\n\t\'2\'\t\ IF THROUGH DATAFILE\n ? ") ; 
scanf ("*/,c" ,&rl ) ; 
scanf("5(c" , Mummy) ; 

printf("\nWAIT !\nWRITING FILE \"datal . in\"\ii") ; 
if (rl==’20 

printf ("\nENTER THE NAME OF DATA-INPUT FILE : ") ; 

scanf ("'/,s" .datain ) ; 

scanf ("*/,c" .Mummy) ; 

ptdatainl = fopenCdatain, "a+") ; 

/* systemC'cp datal.in d.in"); */ 

> 

/♦original = fopen("original","w") ; */ 

/* if(rl=='20 

f printf (ptdatainl, "\n\nENTER THE VALUES OF THE PARAMETERS IN 
THE STRICT ORDER GIVEN\n"); 

f printf (ptdatainl, "BELOW. DATA CAN BE ENTERED IN A LINE ONE 
AFTER THE 0THER\n"); ' 

f printf (ptdatainl, "SEPARATED BY SPACE. \n"); 

f printf (ptdatainl, "DATA SHOULD BE ENTERED ONLY ABOVE THIS TEXT. 
\n") ; 

f printf (ptdatainl, "MORE THAN ONE DATA SETS CAN BE ENTERED- BEGI 
NING OF A NEW\n") ; 

f printf (ptdatainl, "SET SHOULD BE STARTED FROM A NEW LINE FOR 
CLARITY. \n\n"); 

f printf (ptdatainl, "OUTPUT FILES WILL BE GENERATED CORRESPONDING 
TO EACH\n"); 

f printf (ptdatainl, "SET OF DATA.\n\n"); 
fprintf (ptdatainl,"!. VALUE OF SIGMA\n"); 

f printf (ptdatainl, "2. VALUE OF YLst (dimensionless distance)\n") ; 
fprintf (ptdatainl, "3. h, STEP INCREMENT IN THE INDEPENDENT 
VARIABLE TIME\n"); 

fprintf (ptdatainl, "4. TMAX, MAXIMUM VALUE OF INDEPENDENT VARIAN 
LE UPTO WHICHXnDIFF. EQ. IS TO BE INTEGRATED\n") ; 
f close(ptdatainl) ; 

} */ 

printf ("\n\nTWO TYPE OF OUTPUT FILES WILL BE GENERATED BY THIS 
PROGRAM\n\n\ 

1. HAVING NAME AS d(n).o WHERE 'n' DENOTES THE NO. OF THE FILE 
.\n\ 

THIS FILE CONTAINS VARIATION OF DIMENSIONLESS OUTLET CONC. 
WITH\n\ 

DIMENSIONLESS TIME.\n\ 

2. t(n).l THIS FILES GIVES VALUES OF THE TERMS tl, t2, t3, t4 
AND fn\n\ 

USED IN THE PROGRAM . \n\n" ) ; 
if(rl=='2') 
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printfC'NnA FILE NAMED datal.in HAS BEEN GENERATED BY THIS 
PROGRAM FOR ENTERING\n\ DATA IN THE DESIRED FORMAT. Xu 

\nENTER\n\t\'l\'\tIF WANT TO EXIT THE PROGRAM ANDXnXt 
XtEDIT THE ABOVE DATA FILE FOR ENTERING DATA.\a\n\t\' 2\ 
'XtIFX YOU HAVE ALREADY EDITED THE INPUT DATAFILE ANDXn 
XtXtARE READY FOR RUNNING THE PROGRAMXeXh ? ") ; 
scanf ("'/,c" ,&:r2) ; 
scanf ("54c" ,&diiinmy) ; 
if(r2=='20 

/*system("cp d.in datal.in"); */ 

printfC'XnENTER THE NO. OF THE SETS OF DATA : ") ; 

scanf ("*/,d" jftns) ; 

sccinf ("*/,c" , Mummy) ; 

ptdatain = fopen(datain, "r") ; 

printfC'XnXnWAIT ! \n\nSOLVING DIFFERENTIAL EQUATION !\n\n"); 
f or(i=0 ; i<ns ; i++) 

{ 

n=i; 

f scainf (ptdatainl , "*/,g’/,g*/tg54g" ,Asigmajc,ftylst .M.&tmax) ; 
solvelC) ; 

} 

f close(ptdatainl) ; 

> 

> 

if (r2=='l’) 

system("rm -f d.in"); 
exit(-l); 

> 

if(rl==’l’) 

n=0 ; 
do 

printfC'XnXnENTER SIGMA, DIMENSIONLESS REACTION MODULOUS : ") ; 
scanf ("54f" ,Migmax) ; 
scanf ("54c" , Mummy) ; 


printf ("\n\nENTER YL*, DIMEINSIONLESS DISTANCE : ") ; 
scanf ("54f",*ylst); 
scanf ("54c'', Mummy) ; 

printfC'XnXnENTER h, STEP INCREMENT IN THE INDEPENDENT VARIABLE 


TIME : "); 
scanf ("54f",M) ; 


scanf ("54c", Mummy) ; 
("XnXnENTER THE MAXIMUM 
XnWHICH DIFF. EQ. IS 
scanf ("54g" ,&tmax) ; 
scanf ("54c" , Mummy) ; 


VALUE OF THE INDEPENDENT VARIABLE UPTO 
TO BE INTEGRATED : "); 
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scanf ("Xf" j&tmin) ; 
scanf ("*/,c" , Mummy) ; 

printf("\n\uTHE DEPENDENT VARIABLE : ; 

scanf ("y.f",&xin) ; 
scanf ("%c" , Mummy) ; */ 

solvelO ; 

printf("\uDO YOU WANT TO ENTER ANOTHER SET OF DATA ? (y/n) : "); 

scanf ("'/,c" ,&rl) ; 

scanf ("Xc" , Mummy) ; 

if(rl=='y') ++n; 

ns=n+l; 

y while (rl=='y’); 

} 

} 

void solvelO 

■c 

ptdataoutCn] = fopen(f ilenameCn] .fnd,"w") ; 
pttextout2 [n] = fopen(f ilename[n] .fnt2,"w") ; 

/♦fprintf (ptdataout [n] , "sigmax=yf\nylst=y,f\n\n" ,sigmax,ylst) ;*/ 
pttextoutl [n] = fopenCfilenameCn] .fntl,"w") ; 

fprintf (pttextoutl [n] , "♦♦*************************=t‘********\n") ; 
fprintf (pttextoutl [n] , "sigmax=y,f \nylst=y,f\n\n" , sigmax,ylst) ; 
fprintf (pttextoutl [n] , "\ttl\t\tt2\t\tt3\t\tt4\t fuXn") ; 

fprintf (pttextoutl [n] , "*********Nt*********>it**************\ii\n") ; 
tsti=tmin; 
x=xin; 

beta=bel(x) ; 
psi=exp(-beta * ylst) ; 

fprintf (ptdataout [n] , "Xf \ty,f \n" , tst i ,psi) ; 
fprintf (pttextout2 [n] , "’/.f \t%f \n" ,tsti ,x) ; 
j = 1; 
do 

t st i=tmin+h* j ++ ; 
xi=x; 

rkO = h ♦ fnl(tsti,xi) ; 
if ((xi+0.5*rk0)>=0.999) 
breaJc; 

rkl = h * fnl(tsti+0.5*h,xi+0.5*rk0) ; 
if((xi+0.5*rkl)>=0.999) 
break; 

rk2 = h * fnl(tsti+0.5*h,xi+0.5*rkl) ; 
if ( (xi+rk2) >=0 . 999) 
break; 

rk3 = h * fnl(tsti+h,xi+rk2) ; 

X = xi+(rk0+2.*rkl-»-2.*rk2+rk3)/6. ; 
if (xi>=0.999) 
break ; 
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beta=bel(x) ; 
psi=exp(-beta * ylst) ; 
fprintf (pttextout2[n] , "'/.f\t'/,f\n" ,tsti ,x) ; 
fprintf (ptdataoutCn] ,"'/,f\ty,f\]i",tsti,psi) ; 

} while ((tsti<=tm 2 Lx) kk (x<=0.9999)) ; 
f close (ptdataout [n] ) ; 
f close (ptt extout 1 [n] ) ; 

} 

/ ♦♦♦♦*>l' + *******+***************5|t*!*c*****^t***:(c*********!|t******J(<***/ 

float fnl(al,a2) 
float al,a2; 

{ 

float a3; 
float bel(); 
beta = bel(a2); 

a3 = beta * exp(-beta * ylst) + al*0.0; 
return(a3 ) ; 

} 

float bel(a) 
float a; 

{ 

float al; 

sigmajcsql = sigmax * sigmeix; 
tl = pow((l.-a) ,-l./3.) ; 
t2 = log(sigmax)/1.08; 
t3 = t2*t2; 
t4 = exp(-0.9*t3) ; 

al = 1 ./(I ./3 .*tl*tl+2 .♦sigmaxsql*(tl-l . ) + (0.21-0 .62*a)*(l .+ 
sigmajcsql)*t4) ; 

fprintf (pttextout l[n] , "\n '/f \t'/f\t'/Cf Xt'yif XtXf " ,tl ,t2 ,t3 ,t4, al 

); 

return(al) ; 

} 

/♦♦♦♦♦♦*♦♦*♦♦♦**♦♦*♦*♦♦♦*♦*♦*♦♦***♦*****************************/ 
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Sample Simulations 


Table B.l gives the simulations for different assumed values of dimensionless 
reactin modulous, cr keeping dimensionless distance constant (?/£ = 1) 

Table B.2 gives the simulations for different assumed values of dimensionless 
distance, keeping dimensionless reaction modulous constant (cr = 1) 
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Table B.l; Simulations for different assumed values of V’ (yi, - 1) 
Dimensionless outlet cone, for different values of V 


0.0 

0.064 

0.1 

0.065 

0.2 

0.067 

0.3 

0.068 

0.4 

0.070 

0.5 

0.072 

0.6 

0.074 

0.7 

0.076 

0.8 

0.078 

0.9 

0.081 

1.0 

0.083 

1.1 

0.086 

1.2 

i 0.089 

1.3 

0.093 

1.4 

0.097 

1.5 

0.101 

1.6 

0.105 

1.7 

0.111 

1.8 

0.116 

1.9 

0.123 

2.0 

0.130 

2.1 

0.138 

2.2 

0.147 

2.3 

0.157 


0.169 

0.182 

0.198 

0.216 

0.238 

0.263 

0.293 

0.330 

0.373 

0.426 

0.489 

0.56 

0.651 


0.114 
0.114 
0.113 
0.113 
0.113 
0.113 
0.114 
0.115 
0.116 
0.117 
0.119 
0.121 
0.123 
0.126 
0.130 
0.134 
0.139 
0.145 
0.152 
0.160 
0.170 
0.182 
0.195 
0.212 
0.232 
0.257 
0.286 
0.322 
0.365 
0.417 
0.479 
0.550 
0.629 
0.715 
0.802 
0.884 
I 0.952 


0.301 0.354 
0.297 0.356 
0.294 0.359 
0.292 0.364 
0.291 0.370 
0.292 0.378 
0.294 0.388 
0.298 0.400 
0.303 0.414 
0.311 0.430 
0.321 0.447 
0.333 0.466 
0.348 0.487 
0.366 0.508 
0.387 0.531 
0.410 0.554 
0.436 0.578 
0.465 0.602 
0.496 0.626 
0.528 0.649 
0.562 0.671 
0.596 0.693 
0.629 0.714 
0.663 0.734 
0.695 0.753 
0.726 I 0.771 


0.399 0.464 
0.411 0.502 
0.424 0.537 
0.439 0.569 
0.455 0.597 
0.472 0.623 
0.491 0.646 
0.510 0.667 
0.529 0.686 
0.550 0.704 
0.570 0.720 
0.590 0.734 
0.610 0.748 
0.629 0.760 
0.648 0.771 
0.666 0.782 
0.684 0.792 
0.701 0.801 
0.717 0.809 
0.733 0.817 
0.747 0.824 
0.761 0.831 
0.774 0.838 
0.787 0.844 
0.798 0.850 


6.0 10.0 

0.363 

U.Z34 

0.571 

0.698 

0.669 

0.791 

0.725 

0.833 

0.762 

0.858 

0.789 

0.874 

0.809 

0.886 

0.824 

0.896 

0.837 

0.903 

0.848 

0.909 

0.857 

0.915 

0.865 

0.919 

0.872 

0.923 

0.878 

0.926 

0.883 

0.930 

0.888 

0.932 

0.892 

0.935 

0.896 

0.937 

0.900 

0.939 


0.756 0.788 

0.810 

0.855 

0.783 0.804 

0.820 

0.860 

0.810 0.819 

0.830 

0.865 

0.834 0.833 

0.839 

0.870 

0.857 0.847 

0.848 

0.874 

0.879 0.860 

0.857 

0.878 

0.900 0.872 

0.865 

0.882 

0.919 0.883 

0.872 

0.886 

0.937 0.894 

0.880 

0.890 

0.954 0.904 

0.887 

0.893 

i 0.971 0.914 

0.893 

0.897 


0.903 0.941 
0.906 0.943 
0.909 0.944 
0.912 0.946 
0.914 0.947 
0.917 0.948 
0.919 0.950 
0.921 0.951 
0.923 0.952 
0.925 0.953 
0.927 0.954 
0.928 0.955 
0.930 0.956 
0.931 0.956 
0.933 0.957 
0.934 0.958 
0.935 0.959 
0.937 0.959 
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Table B.2: Simulations for different assumed values of ((J - 1) 




Appendix C 


Data Used for Simulation 


Source : All the data given here have been taken from, a Ph.D. thesis (Bhaskar, 
1994). 


C.l Data for active sodium carbonate 


Length of the packed bed ( L ) 

Radius of the pellet (Rg) 

Molar density of the solid reactant { Pp ) 
Grain shape factor { Fg ) 

Void fraction of the bed ( ) 

Space velocity { U ) 

Fraction of reactive solids by weight ( w ) 


= 5.0 cm 
= 0.635 cm 
= 0.0239 g-mole/c77z^ 
= 3 
= 0.42 
= 900 hr-'^ 

_ _i_ 

“■ 40 
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Table C.l: Data related to the active Na2C0z pellet 


Decompn 
Temp., T 

K 

Pore volume, Vg 

( cmVp ) 

Porosity, Cp 

Grain surface 
Area, Ag 

( ) 

Fractional 
vol. of 
sorbent in 
the bed 

423 

0.49 

0.55 

3.14 

0.057 

473 

0.47 

0.54 

2.70 

0.056 

523 

0.45 

0.53 

2.50 

0.055 

573 

0.40 

0.50 

2.20 

0.051 

623 

0.20 

0.34 

0.91 

0.039 


V ariation of rate constant with temperature 


k = A e RT 
where, 


(C.l) 


A' = Frequency factor = 0.022 cm/s 

ka = Activation energy = 4224 cal/g-mole 


V ariation of effective diffusivity with temperature 

D, = 1.76 X 10"“ (C-2) 

C.2 Data for calcined magnesium carbonate 


Length of the packed bed ( L ) 

Radius of the pellet (Ro) 

Molar density of the solid reactant ( Pp ) 
Grain shape factor ( Fg ) 

Void fraction of the bed ( ) 

Space velocity ( U ) 

Fraction of reactive solids by weight ( oj ) 


= 5 cm 
= 0.075 cm 
= 0.0305 g-mole/cm^ 
= 3 
= 0.42 
= 600 hr-^ 

_ _L 
40 


’^TRAL LieRARY 

i 1. T.. KANPUft 






wm 
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Table C.2: Data related to the calcined MgCOz pellet 
(Calcination temp. = 750 °C) 


Pore volume, Vg 

( cm^g ) 

Porosity, Cp 

Grain surface 

Area, Ag 

( ) 

Fractional 
vol. of 
sorbent in 
the 

packed bed 

1.25 

0.82 

11.30 

0.141 


Variation of rate constant with temperature 


k = A e R'T 
where, 


(C.3) 


A' = Frequency factor = 7.12x10 ^ cm/s 

ka = Activation energy = 8213 cal/g-mole 


Variation of effective dif fusivity with temperature 


De = 1 X 10-"® r"®" 


(C.4) 
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Simulated and Experimental 
Results 


Source : All the experimental data given here have been taken from a Ph.D. 
thesis (Bhaskar, 1994). 
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D.l For sulphation of active sodium carbonate 


Sulphation at different temperatures and an inlet cone, of 6000 ppm 


Table D.l: T = 150 ^ (7 


Dimensionless 
Time t* 

Dimensionless Outlet Cone. 

Experimental 

Predicted 
(T=UA,yl = 0.19 

0.006 

0.157 

0.384 

0.011 

0.462 

0.512 

0.022 

0.678 

0.640 

0.033 

0.791 

0.706 

0.045 

0.820 

0.747 

0.056 

0.873 

0.778 

0.067 

0.888 

0.798 

0.083 

0.876 

0.820 

0.100 

0.812 

0.837 

0.117 

0.885 

0.850 

0.134 

0.885 

0.861 

0.150 

0.888 

0.869 

0.167 

0.906 

0.877 

0.184 

0.864 

0.883 

0.200 

0.892 

0.888 
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Table D.2: T = 200 ° C 


Dimensionless 

Time 

Dimensionless Outlet Cone. 

Experimental 

Predicted 
a = 15 A, yl = 0.28 

0.009 

0.131 

0.221 

0.017 

0.426 

0.369 

0.034 

0.626 

0.549 

0.051 

0.683 

0.641 

0.068 

0.747 

0.696 

0.085 

0.753 

0.733 

0.102 

0.730 

0.759 

0.127 

0.753 

0.788 

0.153 

0.761 

0.809 

0.178 

0.822 

0.825 

0.204 

0.788 

0.838 

0.229 

0.788 

0.848 

0.255 

0.805 

0.857 

0.280 

0.813 

0.864 

0.306 

0.826 

0.871 
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Table D.3: T = 250 ° C 


Dimensionless 

Time 

Dimensionless Outlet Cone. 

Experimental 

Predicted 
cr= 15.6, = 0.42 

0.013 

0.103 

0.058 

0.025 

0.153 

0.123 

0.051 

0.361 

0.331 

0.076 

0.465 

0.481 

0.101 

0.572 

0.573 

0.126 

0.634 

0.633 

0.152 

0.685 

0.676 

0.189 

0.712 

0.719 

0.227 

0.685 

0.750 

0.265 

0.715 

0.774 

0.303 

0.754 

0.792 

0.341 

0.768 

0.806 

0.379 

0.805 

0;818 

0.417 

0.826 

0.828 

0.455 

0.838 

0.837 
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Table D.4: T = 300 ° (7 


Dimensionless 

Time 

Dimensionless Outlet Cone. 

Experimental 

Predicted 
(T = 16.4, t/£ = 0.59 

0.018 

0.040 

0.003 

0.036 

0.122 

0.004 

0.071 

0.282 

0.008 

0.107 

0.458' 

0.018 

0.142 

0.407 

0.052 

0.178 

0.504 

0.144 

0.213 

0.496 

0.276 

0.267 

0.606 

0.445 

0.320 

0.648 

0.548 

0.374 

0.680 

0.616 

0.427 

0.689 

0.662 

0.480 

0.712 

0.696 

0.534 

0.757 

0.723 

0.587 

0.803 

0.743 

0.640 

0.821 

0.760 
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Table D.5: T = 350 ^ C 


Dimensionless 

Time 

Dimensionless Outlet Cone. 

Experimental 

Predicted 
a = 17.4,j/£ = 0.66 

0.020 

0.241 

0.0005 

0.040 

0.415 

0.0005 

0.079 

0.612 

0.0006 

0.119 

0.687 

0.0007 

0.159 

0.747 

0.0008 

0.198 

0.747 

0.0010 

0.238 

0.756 

0.0012 

0.297 

0.775 

0.0018 

0.357 

0.834 

0.0030 

0.416 

0.839 

0.0059 

0.476 

0.820 

0.0151 

0.535 

0.855 

0.0513 

0.595 

0.865 

0.1693 

0.654 

0.912 

0.3293 

0.714 

0.934 

0.4548 
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Sulphation of active sodium carbonate at different 
concentrations ( T = 200 °C ) : 


Table D.6: Sulphation at different concentrations 


Time 

(min) 

Dimensionless Outlet Cone. 

1500 

ppm 

3500 

ppm 

6000 

ppm 

18000 ppm 

Expl. 

Pred. 

ri 

Expl. 

Pred. 

t*2 

Expl. 

Pred. 

r3 

Expl. 

Pred. 

t*i 

5 

0.015 

0.104 

0.103 

0.184 

0.123 

0.283 

0.612 

0.579 

10 

0.092 

0.152 

0.426 

0.317 

0.453 

0.467 

0.863 

0.721 

20 

0.233 

0.286 

0.626 

0.512 

0.740 

0.643 

0.910 

0.813 

30 

0.358 

0.388 

0.683 

0.614 

0.880 

0.719 

0.958 

0.852 

40 

0.417 

0.467 

0.747 

0.675 

0.953 

0.764 

0.952 

0.874 

50 

0.417 

0.527 

0.753 

0.715 

0.941 

0.792 

0.943 

0.888 

60 

0.512 

0.579 

0.730 

0.744 

0.956 

0.813 

0.964 

0.899 

75 

0.501 

0.631 

0.753 

0.775 

0.976 

0.836 

0.960 

0.911 

90 

0.565 

0.666 

0.761 

0.798 

0.945 

0.852 

0.976 

0.920 

105 

0.611 

0.696 

0.822 

0.815 

0.964 

0.864 

0.940 

0.926 

120 

0.678 

0.719 

0.788 

0.829 

0.942 

0.874 

0.990 

0.931 

135 

0.717 

0.738 

0.788 

0.840 

0.969 

0.882 

0.965 

0.936 

150 

0.717 

0.754 

0.805 

0.849 

0.901 

0.888 

0.971 

0.939 

165 

0.724 

0.768 

.813 

0.857 

0.957 

0.894 

0.982 

0.942 

180 

0.748 

0.779 

0.826 

0.864 

0.961 

0.900 

0.974 

0.945 


For predicted data, 

(T = 15.4 , yl = 0.27 
r 1 = 1.212 X 10-5 t 
t* 2 = 2.827 X 10-® t 
r 3 = 4.846 X 10'® t 
t* 4 = 14.539 X 10-® t 
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D.2 For sulphation of calcined magnesium carbon- 
ate 

Sulphation at different temperatures [Cone. = 9000 ppm ). 
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Table D.8: T = 750 ^(7 


Dimensionless 

Time 

Dimensionless Outlet Cone. 

Experimental 

Predicted 
<T = 7.5,yl = 0.2i 

0.0056 

0.233 

0.238 

0.011 

0.250 

0.269 

0.022 

0.427 

0.333 

0.034 

0.486 

0.394 

0.051 

0.545 

0.464 

0.067 

0.567 

0.517 

0.084 

0.578 

0.561 

0.101 

0.629 

0.596 

0.118 

0.638 

0.625 

0.135 

0.616 

0.650 

0.168 

0.722 

0.687 

0.202 

0.773 

0.716 


Table D.9: T = 850 


Dimensionless 

Time 

Dimensionless Outlet Cone. 

Experimental 

Predicted 
cr = 0.2, yl = 0.49 

0.012 

0.246 

0.194 

0.024 

0.336 

0.209 

0.047 

0.471 

0.239 

0.071 

0.537 

0.270 

0.106 

0.604 

0.314 

0.141 

0.604 

0.355 

0.212 

0.648 

0.428 

0.247 1 

0.655 ! 

0.459 

0.283 

0.688 

0.487 

0.318 

0.715 

0.512 

0.353 

0.704 

0.534 

0.389 

0.808 

0.555 

1 0.424 

0.814 

0.573 
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Sulpkation at different concentrations (T = 750° C ) 


Table D.IO: Sulphation at different concentrations ( T = 750° C ) 


Time 

(min) 

Dimensionless Outlet Cone. 

2000 ppm 

9000 ppm 

22000 ppm 

Expl. 

Predicted 

ri 

Expl. 

Predicted 

t*2 

Expl. 

Predicted 

t*3 

5 

0.054 

0.217 

0.233 

0.243 

0.412 

0.361 

10 

0.129 

0.227 

0.250 

0.288 

0.583 

0.460 

20 

0.136 

0.245 

0.427 

0.370 

0.728 

0.583 

30 

0.168 

0.265 

0.486 

0.438 

0.796 

0.654 

45 

0.198 

0.291 

0.545 

0.514 

0.781 

0.718 

60 

0.249 

0.318 

0.567 

0.571 

0.837 

0.757 

75 

0.275 

0.343 

0.578 

0.614 

0.813 

0.785 

90 

0.295 

0.366 

0.629 

0.647 

0.826 

0.805 

105 

0.371 

0.389 

0.638 

0.673 

0.852 

0.820 

120 

0.371 

0.410 

0.616 

0.695 

0.881 

0.833 

135 

0.420 

0.429 



0.714 

0.891 

0. 843 

150 

0.450 

0.448 

0.722 

0.729 

0.899 

0.852 

165 

0.491 

0.465 



0.742 

0.906 

0.859 

180 

0.491 

0.481 

0.773 

0.754 

0.917 

0.866 


For predicted data, 

cr = 7.5, yl = 0.24 
t* 1 = 0.694 X 10-® t 
r 2 = 3.125 X 10-® t 
3 = 7.638 X 10-® t 




